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ABSTRACT 



We investigate the number and type of pulsars that will be discovered with the low-frequency radio telescope LOFAR. We consider 
different search strategies for the Galaxy, for globular clusters and for other galaxies. We show that a 25-day all-sky Galactic survey 
can find approximately 900 new pulsars, probing the local pulsar population to a deep luminosity limit. For targets of smaller angular 
size such as globular clusters and galaxies many LOFAR stations can be combined coherently, to make use of the full sensitivity. 
Searches of nearby northern-sky globular clusters can find new low luminosity millisecond pulsars. Giant pulses from Crab-like 
extragalactic pulsars can be detected out to over a Mpc. 

Key words, pulsars: general — telescopes — surveys 



1. Introduction 

Since the discov ery of the first four p ulsars with the Cambridge 
radio telescope dHewish et al.l I1968I) . an ongoing evolution of 
telescope systems has doubled the number of known radio pul- 
sars roughly every 4 years: An evolution from a large fiat re- 
ceiver with a fixed beam on the sky (the orig inal Cambridge 
radio telescope) to focusing dishes (Arecibo - Hul se & Tavloii 
Il975l) . often steerable (Green Ba nk Telescope), on both hemi- 
spheres (the Parkes telescope - Manchest er et alj [20 01 ). with 
large bandwidths and multiple simultaneously usable receivers 
for wider fields of view (Parkes, Arecibo). The types of pulsars 
discovered have changed accordingly, from slow, bright, single 
and nearby pulsars (the original four) to fast (young and millisec- 
ond pulsars), far-away (globular clusters) or dim pulsars, some 
of which are in binaries. 

The next step in radio telescope evolution will be the use 
of large numbers of low-cost receivers that are combined in- 
terferometrica lly. These t elescopes, the All en Telescope Array 
jBoweJ l2007b . LOFAR dRottgerind 120031) . MeerKA T dJonasj 
120071) . ASKAP dJohnston et all 120081) and the SKA dKramej 
2003), create new possibilities for pulsar research. 

In this paper, we investigate the prospects of finding radio 
pulsars with LOFAR, the LOw Frequency ARray. We outline 
and compare strategies for targeting normal and millisecond pul- 
sars (MSPs), both in the disk and globular clusters of our Galaxy, 
and in other galaxies. 



2. LOFAR - The LOw Frequency ARray 

With the first test station operational and the first pulsars de- 
tected LOFAR is on track to start operation in 2010. We have 
evaluated and simulated the LOFAR reference configuration for 
pulsar searches, and will describe that configuration in some de- 
tail below. 

Using two different types of dipoles, LOFAR can observe 
in a low and a high band that range from 30-80 MHz and 110- 
240 MHz respectively. The sensitivity using the high-band an- 
tenna (HBA) is several times that of the low-band antenna (LBA) 



although their survey speeds are similar due to the larger LBA 
field of view. The low band is expected to be a n exciting new 
windo w for exploring radio pulsar behavior (cf. IStappers et al.l 
2007, for an overview of the possibilities for emission physics 
and interstellar medium studies), but the impact on a pulsar sur- 
vey of some of the smearing effects further discussed is so strong 
in the LBA band, that we will only discuss the HBA half of 
LOFAR in this paper. 

The basic collecting elements are the individual dual- 
polarization dipoles; each 4x4 set of these dipoles is combined 
in an analog beamformer and forms an antenna 'tile'. Tiles are 
grouped together in stations; stations farther from the array cen- 
ter are larger than inner stations. In the core, HBA stations are 
grouped in pairs. The innermost 12 stations are 24 tiles each and 
are packed tightly in a 'superstation'. Spread over the 2-km core 
there are 24 more HBA stations of 24 tiles each (making for 36 
core stations in total). These core stations are 32 meters in diam- 
eter, but are tapered to an about 30 meter effective diameter to 
reduce sidelobes. Next there are 18 Dutch 'remote' stations that 
are outside the core and consist of 48 tiles each, while the ~8 
international stations that are spread over Europe use 96 HBA 
tiles. 

At each station the tiles are combined to form up to 8 in- 
dependently steerable 'station beams'. With the cumulative data 
rate out of the station being the limiting factor, the product of the 
number of beams times their bandwidth cannot exceed 32 MHz 
(potentially 48 MHz), e.g.: a station beam set up can range from 
having a single full-bandwidth station beam to 8 independent 
beams of 4 MHz each. Station beams are subdivided in 195 kHz 
channels, and sent to the central processor (CEP) supercomputer 
for correlation, addition and/or different types of beam forming. 
As illustrated in Fig. Q] CEP can further combine station beams 
to form ~ 128 16-bit full-polarisation tied-array beams for the 
superstation, the core and/or the entire array. That number of 
formed beams is limited by the maximum CEP output rate, cur- 
rently 50 Gigabits/second. 

The high-band antennas operate in a 11 0-240 MHz fre- 
quency range and are spaced to optimize sensitivity for the low 
end. They are maximally sensitive toward the zenith. The fall- 
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off of effective area A e ff with zenith angle and observing fre- 
quency has been w ell characterized as shown in Fig. [2] (after 
Wilhelmsson 2007) and we will also describe it qualitatively 
here: at its maximum at 120 MHz, the effective area toward 
zenith is 28 m 2 per antenna tile, halving at a zenith angle of 
~50°. Toward higher frequencies both quantities decrease un- 
til at 240 MHz the maximum effective area is 8 m 2 per antenna 
tile, dropping to half that at ~30° away from zenith. For the 
N s tation = 36 stations of the compact core at 120 MHz, with their 
Nfi/es = 24 tiles each, this effective area translates into a theoreti- 
cal maximum core gain G max ,core of 
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28 m 2 N. 



stations tiles 



2k 



x 10 



-26 



Ws 
Jy m 2 



8.8K/Jy (1) 



where k is Boltzmann's constant in W s m 2 . Finally, the 
noise temperature T an , ema of the HBA is expected to be 14 K at 
120 M Hz and around 180K in the upper half of the band dGunstl 
120071) . 
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Fig. 2. Simulated effective area per antenna tile versus zenith 
angle, for 4 different observing frequencies. From Wilhelmsson 
(200%. 



3. Pulsar searches 

3.1. All-sky surveys 

Although a pulsar survey generally addresses several science 
questions in parallel, it can be optimized for a specific goal: by 
using short integration times potential acceleration smearing is 
kept to a minimum to optimize for finding millisecond pulsars; 
when aiming to maximise the total number of new pulsars found, 
a survey should generally focus on the Galactic plane; while for 
a representative understanding of the local population an all-sky 
survey that is equally sensitive in all directions is optimal. Here 
we will focus on such an all-sky survey. 

Below we first investigate different beam forming scenarios. 
These we compare assuming use of a single station beam, at 
full bandwidth, at the lowest observing frequency of 120 MHz. 
We investigate frequency and band width dependencies in more 
detail in subsection |3.1.2| 
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Fig. 1. Illustration of station, superstation and full-array beam 
patterns at 120 MHz. The zenith is at 0°. The dashed envelope 
is the tile response as shown in Fig. [2] The gains of the different 
beams are here scaled to be the same. Three stations beams are 
shown, up to eight can be formed at a given time. 



3.1 .1 . Beam forming 

To first order an all-sky pulsar survey is optimized to reach best 
sensitivity per given overall time. For that one should minimise 
system noise while maximising collecting area, bandwidth, and 
integration time per pointing. 

For a given survey time and sky area, maximising the time 
per pointing is equivalent to maximising the instantaneous field 
of view. This can be done by forming multiple simultaneous 
beams; while for sparse telescopes like LOFAR, where receivers 
are spread out over a large area, one can also add stations inco- 
herently instead of coherently. This increases the field of view 
at the c ost of increasin g uncorrelated noise over correlated sig- 
nal (cf. lBackerlll999l) . In the incoherent case the beam on the 
sky will be larger, but in the coherent case the system is more 
sensitive. 

From the above we can construct the followin g survey figure 
of merit FoM for a given observing frequency (see Cordes 2002; 
Smi ts et al.ll2009l for related definitions): 



FoM = 10" 



T 

* svs 



beams ^ B 
1 \N stations 



1/2 



(2) 



where the factor 10 3 scales the reference FoM to be 1.0; 
A/T sys is the ratio of effective area and system noise equiva- 
lent temperature, averaged over the field of view, a ratio indicat- 
ing telescope gain. Nbeams is the number of simultaneous beams 
formed and O denotes the field of view as derived from the beam 
size. Here and below we estimate all beam sizes by their full 
width at half maximum (FWHM) as 1.22^deg with A the 
observing wavelength, and D the diameter of the element over 
which the beam is formed (e.g.: a station, the core, or the full 
array). Then Q = n( FW ^ M ) 2 deg 2 . B is the observing bandwidth; 
the right-hand denominator is 1 for coherent addition, or N stat i ons 
when that number of stations is added incoherently. This FoM 
is inversely proportional to the minimum detectable flux S m ,„. 
When comparing between setups, the ratio of survey speeds - 
the time needed to reach a given sensitivity over the whole sky 
- is FoM 2 . In Table Q] we list the FoM for three survey set ups 
discussed below. 

To outline the various trade-offs involved with these survey 
setups, we first compare scenarios for coherent versus incoherent 
addition for the LOFAR core only, and not the total array - given 
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Type 


A/T n . s (nrVK) 


beams 


Si (deg^) 


B (MHz) 


1 l^slations 


FoM 


Full Incoherent 


260 


l 


27 


32 


54 


1.0 


Core Coherent 


170 


128 


0.006 


32 


1 


0.8 


Superstation Coherent 


40 


128 


0.23 


32 


1 


1.2 



Table 1. Figure of merit FoM per Eq.[2]for the three different survey set ups as described in the text: Full Incoherent, Core Coherent 
and Superstation Coherent. All parameters derived using the lowest observing frequency, 120 MHz. In the first column, A/T sys , we 
use 54, 36 and 12 stations respectively, of 24 tiles at 28 m 2 , and 140 K, while for the Superstation Coherent case we also multiply 
by 0.70 to correct for the lower sensitivity at the station beam edge. 



the sparseness of the distribution of remote stations, coherent 
addition of the entire array would result in an extremely limited 
field of view. 

In the Core Coherent scenario all 36 stations in the 2 km di- 
ameter D core compact core are combined coherently. In this case 
Q = 0.0060 deg 2 and the gain is as per Eq. [TJ For comparison, 
each of the individual roughly D sta f wn = 30-meter core stations 
forms a beam on the sky of 27 deg 2 . If added incoherently, the re- 
sulting beam is as large as for an individual core station (cf. Fig. 
Q]) but at a factor y/N sla ,i „ s = 6 decreased sensitivity compared 
to the coherent case. The factor {-§^-) 2 = (^f^) 2 =4.4 x 10 3 
increase in beam area allows for longer integrations by the same 
factor. For synthesis telescopes that form a single coherent 'tied- 
array' beam, like the Westerbork Synthesis Radio Telescope or 
the Very Large Array, such a factor would well describe the 
trade-off; but as the LOFAR CEP can output N beams = ~ 128 co- 
herently added core beams simultaneously versus only a single 
incoherently added core beam, the factor F, difference in inte- 
gration times is F, = (J^-f/N heams = (») 2 /128 = 35. Thus, 

in this comparison the sensitivity decrease by y/N stat ions — V36 
of adding stations incoherently and the sensitivity increase by 
Vf7= V35 provided by longer integrations practically cancel 
out: the FoM for coherent or incoherent addition of just the core 
is very similar. 

Although for the coherent sum the addition of remote sta- 
tions significantly reduces the field of view, for the incoherent 
case the entire array can be added without field-of-view penalty. 
Physically the remote stations are twice the size of the core sta- 
tions but they can be tapered to the size of the core stations to 
produce a similar field of view. Adding these 18 remote stations 
in the above equation produces a FoM for the Full Incoherent 
scenario of 1.0; that is (F t N s , a ,i ons j u iiy /2 /N sla ,ioM,co r e = 05 X 
54y/ 2 /36 = 1.2 times higher than FoM for the Coherent Core 
case which is 0.8 (Table [TJ. Yet (only) for the Coherent Core 
survey producing more simultaneous beams increases the FoM; 
if 50% more beams could be formed these two scenarios would 
be equally sensitive. 

Incoherent addition may offer slightly better sensitivity be- 
cause it allows for long integration times. If however these long 
integration times are a significant part of a potential pulsar binary 
orbit (cf. 1-hr integrations in our reference incoherent survey, 
defin ed below, versus th e 2.4-hr binary period of PSR J0737- 
3039 Bur gav et alj 12003). the apparent change in pulse period 
due to the acceleration through the orbit decreases the effec- 
tiveness of periodicity searches. At a computational cost differ- 
ent types o f acceleration searches can mitigate this decreas e in 
sensitivity dJohnston & Kulkarnil[l99l iRansom et ai1l2003l) . A 
coherent survey is less prone to this acceleration smearing: to 
reach the same S m m as the incoherent reference survey with 1-hr 
pointings, the coherent addition survey integration time is down 
to 3 minutes. From an efficiency point of view short pointings 
are also more robust to system errors such as data glitches, and 



to impulsive radio interference. In contrast to searching a sin- 
gle beam for an incoherent-addition survey, handling ~128 si- 
multaneous data streams, each at high time resolution and total- 
ing 50 Gigabit s/second is a sig nificant computat ional and data- 
handling task dBackerll 19991; cf. lSmits et ai1l2009l for the SKA). 

In a different coherent survey, Superstation Coherent, 128 
tied-array beams from just the superstation (cf. section |2| tile 
out the complete station beam. Compared to coherently adding 
the entire core, the superstation beams have larger field of view 
but are less sensitive. In a grid of 128 of these wider supersta- 
tion beams many will be close to the less sensitive edges of the 
station beam, causing for some decrease in sensitivity compared 
to the narrow Coherent Core beams. For A/T svs in Table [TJ we 
hence use the average sensitivity out to FWHM for an assumed 
Gaussian station beam, which is 70% of the sensitivity for the 
center of the beam. Yet the relatively large field of view caused 
by the dense packing of the superstation offsets the limited sen- 
sitivity of the lower number of stations now used, producing a 
FoM of 1 .2, compared to 1 .0 for the incoherent survey discussed 
above. Given identical overall survey duration, the integration 
times for this coherent superstation and the incoherent full ar- 
ray are similar, so there is no reduction in acceleration smearing. 
Inherent to coherent-addition surveys, the computational task of 
searching through 128 beams simultaneously and in real-time 
remains. 

Beam forming scenarios have an impact on how accurately 
sources can be located; this we will address in more detail in 
section|4]where we discuss follow up. 

3.1.2. Observing frequency 

Determining the optimum observing frequency involves map- 
ping out the trade-off between beam size, source and back- 
ground brightness, and pulse smearing from dispersion and scat- 
tering. Some telescope and pulsar behavior favors low-frequency 
observing: when moving up from the lowest frequencies the 
beam size decreases as v~ 2 which limits integration time per 
pointing in a fixed-time all-sky survey; the effective collecting 
area falls off with frequency; and the intrinsic p ulsar brightness 
decre ases with frequency as v~ L5 on average (Malofeev et al. 
2000). Other circumstances are more favorable at higher fre- 
quencies: the backg round sky noise decreases steeply as v -2 6 
dLawson et al.lll987l) and smearing from interstellar dispersion 
and scattering, which makes the pulsar periodicity harder to de- 
tect, is less. Much of the effect of dispersion can be removed 
by se arching over many (~10 4 ) incoherently dedispersed trial 
DMs (Lorimer& Kramer 2005). No such removal process ex- 
ists for scatter broadening, where the pulsar emission takes dif- 
ferent paths through the interstellar medium with different travel 
times, resulting in an observed pulse profile that is significantly 
smeared out. Scattering smearing time in creases sharply with 
decreasing observing frequency (as v • , iRamachandran et all 
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Fig. 3. Number of pulsars detected in simulated surveys at dif- 
ferent frequencies and sensitivity estimates, plotted as a ratio 
over the productivity of our reference 25-day all-sky survey of 
1-hr pointings at 140 MHz. Towards lower frequencies the yield 
goes up, mostly caused by the increasing effective area as the 
wavelength approaches the telescope's optimum at twice the 
antenna spacing of 1 .25 m; equally important is the increasing 
beam size which allows for longer integrations for the same to- 
tal survey time. 



clear, and RFI is significantly reduced by the 12-bit digitisation 
and the many-element interferometer design, especially when el- 
ements are combined coherently. As the total bandwidth can fur- 
thermore be quickly split up and spaced to avoid channels con- 
taminated by RFI, we expect that the entire specified 32 MHz 
band (potentially 48 MHz) will be usable, but below we shall 
conservatively use 80% of 32 MHz. The potential 1.5-fold in- 
crease in usable bandwidth from 32 to 48 MHz would increase 
sensitivity by V48/32. As S„„„ depend s on zenith di s tance and 
sky noise it varies per pointing, but after lDewev et al.l d 1985b we 
can estimate a typical value for a 1 -minute pointing towards the 
zenith using a coherently formed core beam, for a pulsar with a 
10% duty cycle: 



= SNR 



= 10 



T 

1 sys 



0.66G ma AN pol 0.8Bt)2 



\P-W) 



1.0 x 10 3 K 



v 0.66x8.8^(2x0.8x32MHzx60s)3 J 
10.3 mJy 



(3) 



For comparison the S m ,-„ for incoherently added pointings 
for our reference survey described below (1-hr integrations, 54 
stations added) comes to 5.8 mJy. 



119971; or alternatively v~ 3 - 9 . iBhat et alj|2004l) and with increas- 
ing distance to the pulsar dTavlor & Cordeslll993l) . 

Although the sky background and scatter broadening in- 
crease towards lower frequency we find in our simulations, fur- 
ther described below, that the total survey productivity, if de- 
fined as the total number of pulsars detected (Fig. [3), is mainly 
determined by the effective collecting area and the beam size. 
These peak towards the lower edge of the band, leading us to 
conclude that the survey is most efficient at lowest frequencies. 
When using a single full-bandwidth 32 MHz station beam this 
means a central frequency of 140 MHz, from here on our refer- 
ence frequency. As can be seen from Eq. [2] bandwidth can be 
freely traded for beams with no impact on the FoM. Splitting 
up the available bandwidth over multiple independently point- 
ing stations beams (keeping beams-bandwidth product equal to 
32 MHz as outlined in Section |2]l and moving each beam to a 
lower, more sensitive observing frequency may therefore further 
increase the overall survey output. This does increase integra- 
tion time per pointing, which may be only beneficial in the Core 
Coherent scenario where short default integration times could 
limit sensitivity to intermittent pulsars. 

3.1.3. Sensitivity 

The minimum detectable flux S mi „ depends on the signal-to- 
noise ratio at which one accepts a signal as real (SNR), the 
system temperature (T J> , J =T ree +T J ^) as it varies over the sky, 
the number of polarisations (N po i), the bandwidth (B), the inte- 
gration time (t), the pulsar period (P) and pulse width (W) and 
the zenith-angle dependent gain G(z). Compared to the theoret- 
ical gain the real-life gain is always less; although in Eq. Q] we 
already use the simulated effective area and take aperture effi- 
ciency into account there, we apply a conservative factor 0.66 to 
estimate the real-life gain for our simulations below. Regarding 
the impact of RFI on the amount of usable bandwidth, testing 
has shown the radio-interference environment to be relatively 



3.1.4. Simulations 

To determine the number and type of pulsars LOFAR can find, 
we have used the above characteristics to model the detec- 
tion of a large number of simulated normal pulsars. For this, 
we have used a population synthesis code that simulates the 
birth, evolution, death and possible dete c tion of radio pulsars 
(for details see lBhattacharva etal]|l992l;lHartmanetal.lll997l: 
Ivan Leeuwen & Ver bunt 2004). We do not simulate millisecond 
pulsars, as a realistic treatment of the survey selection effects 
associated with their binary history would be essential but is 
beyond the scope of this wor k. For our simulated p ulsars we 
draw birth velocities from the iLvne & Lorimerl ( 1 19941) distribu- 
tion with its 45 + 90 km s~' mean v elocity. Initial positions are 
simulated as in lHartman et al.ld 19971) : heights above the Galactic 



Survey 


Year 


N real 


N 


sim 


v (MHz) 


Jodrell 


1972 


51 


20 


± 2 


408 


UMass-Arecibo 


1974 


50 


39 


± 4 


430 


MolongloII 


1978 


224 


227 


± 7 


408 


UMass-NRAO 


1978 


50 


54 


± 6 


400 


ParkesII 


1991-1994 


298 


335 i 


10 


400 


Cambridge 80MHz 


1993-1994 


20 


27 


± 5 


82 


Parkes Multibeam 


1997-2000 


987 


801 i 


41 


1374 


LOFAR 


2010-2012 




1100 ± 


100 


140 



Table 2. Name, year, number of detected pulsars (N rea i), 
number of simulated detected pulsars (N s j m ) and frequency 
observed at (v) for the eight surveys simulated. The er- 
ror on N s j m is the standard variation in the outcomes 
when simulating the model described in the text with dif- 
ferent random number seeds. The N ref| i d ata is taken from 
iDavies et alJ d 19771): iHulse & Taylorl dl975h: IManchester et al.l 
d!978l): iDamashek et alj dl978l): IManchester et al J dl996): 
Shra uner et al.l d 1998b and Manch ester et alj (2005) respectively. 
The LOFAR numbers are for the incoherent 25 -day reference 
design described in the text. 
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plane are randomly selected from an exponential distribution 
with a scale height of 60 pc; for the initial distribution in galac- 
tocentric radius we use their Eq. 4, modified to scale per unit 
area: 

p(R)dR= |L expJ-|-JdR (4) 

where scale length R w is 5 kpc. This modified distribution does 
not peak as strongly toward the Galactic centre; for the local 
population there is no significant change but surveys that sam- 
ple the inner Galaxy are better reproduced. We next calculate 
the 3D orbits through the Galaxy. We evolve magnetic fields and 
periods, determine whether the pulsar is still ab ove the death- 
line. W e estimate the luminosity at 400 MHz (per lHartman et all 
119971 who follow iNaravan & Ostrikerll 19901) . With a spectral in- 
dex drawn from a Gaussia n distribution of width 0.76 around 
-1.5 dMalofeev et alJuOOOl takes luminosity turn-overs into ac- 
count) we scale this luminosity to the observing frequency. 

Using the simulated position w e look up the sky background 
temperature (Hasl am et al.l [1982) and combined with the dis- 
tance we also esti mate the dispersion m easure and the scatter 
broadening from iTavlor & Cordesl (119931) . Given that both sky 
background and scatter broadening time increase toward lower 
observing frequencies we next scale these u sing power laws with 
spectral index -2.6 an d -4.4 respectively dLawson et al.l ll987; 
ITavlor & Cordell993h . 

For each of the 8 pulsar surveys in Table|2]we model the sen- 
sitivity function versus period, dispersion and scattering smear- 
ing, and sky position. We also model the decrease of pulse width 
and hence detectability with period. By comparing the simulated 
and real pulsa r sam ples for the first four of the above surveys 
lHartman et al.l d 1997b determined the most probable underlying 
initial model parameters. 

We use this best model to determine the yield of a future 
survey with LOFAR. As tests, we also model the second 81.5- 
MHz Cambridge survey dShrauner et alj[l998l) to check the va- 
lidity of our extrapolatio ns to lower frequencies and the Parkes 
Multi beam pulsar survey (M anchester et alfcOOUlLorimer et alj 
2006, hereafter PMB) because of its superior statistics. These 
tests we describe in some more detail in the two paragraphs be- 
low. 

For each simulated pulsar, we model if it would be de- 
tected by the Cambridge survey given its sky position, disper- 
sion and scatter smearing, and luminosity. We estimate the min- 
imum detectable flux fro m the S m /„ versus DM and period plot 
in lShrauner et al.l ([T998), their Fig. 3. That figure does not take 
scattering into account, although for these simulations one real- 
istically should. For each simulated pulsar we therefore deter- 
mine the 'effective' (higher) DM-value that would produce the 
same smearing as the actual dispersion and scattering do com- 
bined. We next determine the sensitivity for the period and that 
'effective' DM. The sensitivity thus estimated is valid for the sur- 
vey's median sky background noise of 2 Jy, so we next scale it 
using the actual sky noise for the position of the simulated pulsar. 
Finally we multiply the s ensitivity by the declination-dependent 
power response (Fig la., Shraun er et al.l 1 1 998b for the sky po- 
sition and check if the simulated pulsar is bright enough to be 
detected. With 27 ± 5 simulated pulsars found, our model re- 
produces the actual tally of 20 pulsars quite well, validating the 
luminosity, sky-noise and scattering extrapolations to lower fre- 
quencies that we also use for the LOFAR predictions. 

As a further check on our model, the PMB simulations find 
800 ± 40 pulsars, which compares reasonably well to the actual 
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Fig. 4. Simulated number of pulsars detected in an all-sky sur- 
vey, versus total survey time in days. A survey time of 25 days 
corresponds to 1-hr pointings for incoherently added stations 
beams. 



numb er of 987 non-recycled pulsars detected dManchester et al.l 
120051 ATNF catalog per Jan 1, 2009). 

For LOFAR, we have evaluated surveys for different tele- 
scope parameters, total survey time and minimum sensitivity. 
Our reference survey is a Full Incoherent survey using 54 sta- 
tions, 64-minute pointings, a gain per Eq. [3] of 0.66 times the 
theoretical gain, with 80% of 32MHz bandwidth at a central fre- 
quency of 140 MHz. With a single beam such a survey can scan 
the visible sky in 25 days. In an all-sky survey with this reference 
set up 1 100 + 100 pulsars could be detected (the full curve in Fig. 
H), 900 + 100 of which are new. Our sensitivity estimate B uses 
a less conservative sensitivity of 0.8 times the theoretical gain. 
In the same overall time it finds about 100-200 more pulsars 
that are new; i.e. not detected by older surveys we simulate. We 
do not take into account findi ngs from ongoing surveys such as 
the G BT350 search project dHessels et al.ll2007h lArchibald et al.l 
2009). Here, and below, the errors quoted for simulation runs in- 
dicate variations in output caused by using different initial ran- 
dom seeds. Between simulation runs, the number of new detec- 
tions is a stable fraction of the total number of detected pulsars, 
indicating that in the simulations there is a well-defined part of 
the pulsar population that only LOFAR is sensitive too. In con- 
trast, there are also known pulsars that LOFAR will not be able to 
detect in this reference survey, mainly due to scattering effects. 
We have also evaluated pointings of 1,2,4-256 minutes duration 
(see Fig. 0). Here the duration effects the minimum detectable 
flux; using the figures of merit in Table Q] this incoherent survey 
can be scaled to the other survey types. 

If pulsar luminosity is the limiting factor in the detections, a 
larger time per pointing t decreases the minimum detectable flux 
S m in as ?~i , increasing the distance d out to which pulsars can be 
detected as f*. At the scale of the Galaxy pulsars are located in 
a disk; if one could observe pulsars throughout the Galaxy, the 
number of detectable pulsars is N ~ d 2 ~ t* . Locally, pulsars are 
distrib uted roughly isotropically; then N ~ d 3 ~ t* (cf. ICordesI 
2002). We find that for LOFAR surveys, N is even less dependent 
on t than N ~ n because scatter broadening is the limiting factor, 
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Fig. 5. The 1100 simulated pulsars detected with the reference Galactic LOFAR survey, for 1 -hour pointings. The scattering inter- 
stellar matter (ISM) is shown in gray. Left) projected on the Galactic plane. Right) projected on the plane through the Galactic 
centre and sun, perpendicular to the disk. The Galactic centre is at (0,0,0). LOFAR probes the local population very well, while the 
broa dening from multi-path scattering on the ISM limits what volume a low-frequency survey like LOFAR probes. (ISM modeled 
after lTavlor & Cordesll 1993b . 



not luminosity (cf. Fig. |4] where the V- line is scaled to be equal 
to the reference sensitivity estimate at the 6.25 day point). 

Compared to the Parkes Multibeam survey, the pulsars found 
with LOFAR are different in several ways. The PMB is more 
strongly focused on the Galactic centre. The range of the pop- 
ulation detected by LOFAR is limited by scatter broadening. 
Because of its higher sensitivity, LOFAR detects more low- 
luminosity nearby pulsars (Fig.[5k). A lower limit to p ulsar lumi- 
nosity of 1 mJy kpc 2 at 400 MHz has been suggested ( Lvne et al. 
1998, cf. similar lower limit at 1400 MHz in lLorimer et alj 
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2006); with a spectral index of -1.5, this compares to a lumi- 
nosity at 120 MHz of 6mJykpc 2 . If this lower limit exists the 
reference survey can detect all pulsars that are beamed toward 
us up to 1 kpc. If it does not exist, then this survey will cer- 
tainly illuminate how many lo w-luminosity pulsar s were formed 
nearby (cf. PSR J0240+62 in lHessels et al.ll2007l) : a number of 
importance if one wants to understand the neutron-star birthrate. 

Furthermore, LOFAR finds more pulsars out of the Galactic 
disk (Fig. |5j3) as its high sensitivity is unimpeded by scattering 
in that direction. As pulsars are born in the plane (z=0), this 
observed z- distribution could disclose their much deb ated birth 
velocity (cf. Hartmanl l 199711 Arzoumanian et al.l l2002). 

For the detection of millisecond pulsars (MSPs) scatter 
broadening will be the limiting factor. Although there is much 
variation in scattering measure for different lines of sight, on 
average the scatter broadening at 120 MHz is on the order of 
lms for sources with DMs higher than 30pc/cm 3 (using the 
rough empirical t sc =2 x 10 9 D M 38 relation between DM and 
scattering time by IShitovll 19941) or d istances of about 1 kpc (us- 
ing distance-scattering relation from lTavlor & Cordeslll993l) . As 
many of the currently known Galactic-disk MSPs are within this 
range o f 1 kpc in which scatter broadening does not hinder de- 
tection dManchester et al J 120051) . LOFAR could probe the local 
population to a much lower flux limit. 



20 



a 
z 



4 5 6 
Period (ms) 



10 



Fig. 6. Periods of the 115 MSPs in globular clusters with pe- 
riods les s than 10 ms . There are another 23 pulsars with longer 
periods. (lFreirell2008l) . 



3.2. Galactic globular cluster surveys 

In a fixed-time all-sky survey, one can gain field of view at a 
cost of instantaneous sensitivity by adding the signal of the sta- 
tions incoherently instead of coherently. When the incoherent 
FoM is higher (as in Table [D this trade-off increases the number 
of detectable pulsars because it allows for more time per point- 
ing. In contrast, specific, smaller regions on the sky with higher 
densities of radio pulsars can potentially be better targeted with 
smaller field of view, but significantly higher sensitivity. 

Globular clusters fit the description well; they are compact 
and form regions on the sky with high stellar densities. These 
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high densities also cause globular clusters to contain more bina- 
ries and binary products than are found in the disk. This makes 
globular clusters very good candidates for MSP searches (e.g. 
iLvne et aTlll987l) . To estimate which clusters are most promis- 
ing for a LOFAR search, we evaluated several of their proper- 
ties: location on the sky, dispersion measure (DM), distance (d), 
and the number of radio pulsars potentially present. 

As discussed above for sources in the Galactic plane, scatter 
broadening is a concern for detecting far-away fast millisecond 
pulsars. Most globu lar cluster pulsars have periods of around 
2-5 ms dFreird 120081 see Fig. [6j. Using the above mentioned 
DM versus scattering time relations, at 120 MHz such MSPs 
are detectable up to a DM of about 30^-0 pc/cm 3 . Observing at 
200 MHz extends this limit to 60-80 pc/cm 3 . For the longer pe- 
riod pulsars also present in these clusters, these limits are less of 
a problem, so we will further investigate globular clusters with 
DMs up t o 100 pc/cm 3 b e low . From studies with the Westerbork 
telescope Stappers et al. (200§|) also conclude that the DM ver- 
sus scattering time relation is uncertain to the extent that some 
DM = 100 pc/cm 3 MSPs may be detectable. We rate the candi- 
date clusters by the expected number of detectable pulsars: as- 
suming j ffj; — _ 1 dAndersonl 1 992f) . this number scales as d~ 2 , 
as the declination-dependent telescope gain G(z) (cf. Figure [2]), 
and as the collision number F = p\rji, where p ? and r c are th e 
central density and the core radius, respectively (IVerbuntl l2003). 
The sources that are most promising according to this scaling 
are listed in Table [3] The half-mass radius of each fits within a 
LOFAR coherent core beam (cf. Fig.|7]for highest-ranking can- 
didate M15), which means we can now use the full gain of the 
coherent addition. 

In general, searching globular clusters more deeply than 
done previously requires integrating for several hours and 
analyzing the d ata with multiple b inary-acceleration search 
techniques (see iRansom et alj I2005L and references therein) . 
Assuming an MSP spectral index of —1.7 dKramer et al.l fl998) 
the weakest pulsar in M15 (B2 127+1 1G, 0.13 mJy at 430 MHz, 
from Anderso n! 19921) is 1.1 mJy at 140 MHz; in a 10-hour point- 
ing with the compact core at 140 MHz, the minimum detectable 
flux (eq.[3]at M15's zenith angle of 40°) is 0.5 mJy. Thus far, the 
luminosities of the pulsars in M15 follow | '°^ L = — 1. If this re- 
lation extends to lower luminosities, observing up to a minimum 
flux of 0.5 mJy could yield several new MSPs. 



Name 


DEC 


d 


DM 


r c 


r 


N 




deg 


kpc 


pc/cm 3 


arcmin 






M15 


+ 12 


10.3 


67.3 


0.07 


665 


8 


M92 


+43 


8.2 


29+ 


0.23 


106 




M5 


+02 


7.5 


30.0 


0.42 


79 


5 


M10 


-04 


4.4 


43 + 


0.86 


28 




M13 


+36 


7.7 


30.4 


0.78 


39 


5 


M3 


+28 


10.4 


26.3 


0.55 


66 


4 


M2 


-01 


11.5 


29 + 


0.34 


120 




M14 


-03 


9.3 


76 + 


0.83 


54 




Pal 2 


+31 


27.6 


97 + 


0.24 


209 




M12 


-02 


4.9 


39+ 


0.72 


10 





Table 3. Name, declination (DEC), dist ance (d), dispersion mea - 
sure (DM) observed or + modeled after Tay lor & Cordesl (1 19931) . 
core radius (r c ), collision number (T) and number of detected 
pulsars (AO for the ten highest-ranking candidates f or a L OFAR 
globular cluster survey. Data from lHarrisI d 19961) and iFreird 



(2008). 



3.3. Pulsars in other galaxies 

The next over-dense regions on the sky are galaxies. Their dis- 
tance is a problem, but compared to globular clusters, galaxies 
have several advantages for a LOFAR pulsar survey. They are 
distributed equally over both hemispheres, while most globular 
clusters are invisible from the telescope site in The Netherlands; 
and if visible face-on and located in the part of the sky that is 
pointed away from our Galactic disk, the scatter broadening is 
relatively low. 

In the relatively nearby galaxy M31, a periodicity search 
based on a 10 hour pointing with the compact core at 140 MHz 
could detect all pulsars more luminous than ~200Jykpc 2 
(see eq. |3), which is comparable to the top end of the lu- 
minosity distribution in our own galaxy. Any M31 pulsars 
that emit giant pulses dArgvle & Gowerl 1 19721) could be dis- 
covered more easily through these giant pulses than through 
their periodicity if the flux ratio between giant and normal 
pulses exceeds 10 5 dMcLaughlin & Cordesl 120031) . Two young 
pul sars, the Crab pulsar B05 3 1 +21 and LMC pulsar B0 540- 
69 dStaelin & Reifenstein|[T968l; iJohnston & Romanill2003l) . are 
known to emit such giant pulses and the former was indeed dis- 
covered through them. Compared to the regular periodic emis- 
sion, giant pulses show a steeper sp ectral index (-3.0 to -4.5, 
IVoutdl200U iMcLaughlin & Cordesl 120031) making giant-pulse 
searches especiall y attractive for a low-frequenc y telescope like 
LOFAR. Using the lMcLaughlin & Cordesl d2003l) maximum dis- 
tance estimator for giant-pulse detection in 1-hr pointings 



d = 0.85 Mpc 



5Jy 



> CP 



10 5 Jy 



n 



B \ 



1/4 



10 MHz ) 



(5) 



with S sys = G core /T sys = (0.8 x 8.8 K/Jy)/(1.0 x 10 3 K) = 140 Jy, 
bandwidth B-32 MHz and scaling the giant pulse flux with 
spectral index -3.0, the Crab pulsar would be detectable out to 
-1.5 Mpc. 
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21 h 30 m 10 s 21 h 30 m 00 s 21 h 29 m 50 s 
Right ascension (J2000) 

Fig. 7. Beam setup for a MSP survey of M15, the globular clus- 
ter with the highest success probability in Table [3] The cluster 
half-mass radius falls well within t he core beam. T he crosses 
mark the 8 pulsars currently known dAndersonl ll992). 
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Fig. 8. Seven coherently formed beams from the 2-km LOFAR 
core projected on candidate galaxy M8 1 at their FWHM. 



Detection of a handful of pulsars in each of several nearby 
galaxies would enable comparison of the top end of their lu- 
minosity functions, and potentially show a relation with galaxy 
type. Through pulse scattering and dispersion these pulsars will 
also probe intergalactic matter. 

Like for the globular clusters above we have ranked target 
galaxies by distance, size on the sky and expected number of 
sources (here assumed to scale linearly with mass). We select 
galaxies that are within 5 Mpc from the Earth and that are not 
too strongly scatter broadened given their Galactic latitude gb 
and inclination angle i (requiring \gb\+9Q-i > 30). In Table|4]we 
list the 10 highest-ranked candidates. For most of these highest- 
ranked galaxies the core can be captured with a single compact- 
core beam (cf. Fig. [§), or the whole galaxy by a pattern of 7 
beams. 



Name 


d 

Mpc 


s 

arcmin 


i 

deg 


gb 
deg 


logM 

M Q 


M31 


0.7 


193 


78 


-22 


11.4 


M81 


1.4 


22 


60 


41 


10.7 


M33 


0.7 


56 


56 


-31 


10.1 


M94 


4.3 


12 


33 


76 


10.8 


NGC2403 


4.2 


23 


62 


29 


10.7 


NGC4236 


2.2 


19 


73 


47 


9.8 


NGC4244 


3.1 


15 


90 


77 


10.0 


NGC4395 


3.6 


13 


38 


82 


10.1 


NGC3077 


2.1 


5 


43 


42 


9.1 


UGC7321 


3.8 


5 


90 


81 


9.7 



Table 4. Name, distance (d), size on the sky (s), inclination angle 
(0, Galactic latitude (gb) and logarithm of the mass (log M) for 
the ten highest-ranking candi dates for a L OFAR extragalactic 
pulsar survey. Data taken from lTulfvl dl988h . 



4. Discussion 

4. 1 . Steep spectrum sources 

Generally pulsars are increasingly bright toward lower fre- 
quencies, but this steep spectrum of pulsars typically flattens 
out or turns oyer at a frequency between 100 and 250 MHz 
dMalofeev et alj Il994t) . There is, however, a small fraction 
of pulsars for which no such spectral break is observed 
down to frequencies as low as 50 MHz. Sources with a spec- 
tral index steeper than the spectral index of the sky back- 
ground of -2.6 (e.g. PSR B0943 +10, spectral index -4.0; 
Ramachandran & Deshp andel 1 19941) will be more easily de- 
tectable by a telescope like LOFAR, producing quantitative input 
for radio emission models. 



4.2. Follow up 

With high sensitivity most important and field of view not a con- 
sideration, follow up timing (t he determination of puls ar and po- 
tential binary parameters, cf. lLorimer & Kra mer 2005) will use 
the entire core, and potentially more of the LOFAR array if it can 
be reliably phased up, for beam forming. Given the small 3' field 
of view, such follow up needs very accurate candidate source 
positions. Determining the source position that accurately also 
directly facilitates timing follow-up by eliminating position as a 
variable. 

In some of the above surveys sources are only localised to 
within the 30' superstation or 5°. 8 station beams (cf. Section 
13.1. 11 1. Detections can be more accurately located later by fol- 
lowing up with many, smaller tied-array beams. As the DM 
is now known the total computational burden of searching 
through these many tied-array beams is not as problematic as 
for an entirely coherent survey. For example, we find that about 
50% of incoherent-survey pointings will have sources brighter 
than the reference S m i„. Each of these station-beams pointings 
could subsequently be filled out with 128 coherent supersta- 
tion beams in about l/3rd of the original total survey time 
(50% ( F f;' Mf "" f = 0.36, cf. Eq.|2]and Table Q) to reach the 

r OM Super station 

same 5„„„. This would provide the necessary first conformation 
observation but also an improved position. In a next step 128 
compact-core beams can then similarly tile out the superstation 
beam in which the source is located to provide full localisation, 
and the first timing data. 

Once sources are properly localised, long-term follow up 
timing uses tied-array beams. From our simulations we find that 
in the incoherent reference survey about 80% of sources have 
nearest neighbors at less than one FWHM station beam away 
and hence these pairs of sources can be timed simultaneously. 
18% and 12% of simulated pulsars are in groups of 3 and 4 per 
station field-of-view respectively. In total the follow up needs 
to re-observe about 50% of the survey pointings to time each 
newly found pulsar once. Sub-arraying or trading bandwidth 
for independently steerable station beams could potentially in- 
crease overall follow-up timing efficiency. The timing observa- 
tion should hold at least several hundred pulses for a steady pro- 
file to form, which at an average 1 s-pulsar translates to about 
lOmin integrations. Given the yIN stations = 7.3-fold increase in 
gain between incoherent surveying and tied-array follow-up, a 
timing run can thus produce high SNR detections for all candi- 
dates. 



Joeri van Leeuwen and Ben Stappers: Finding pulsars with LOFAR 



9 



4.3. Fast radio transients 

Given the software-telescope nature of LOFAR, piggybacking 
and simultaneous observing are relatively straightforward to im- 
plement. As the 'telescope time' concept at LOFAR includes 
the availability of central signal processing at CEP, observation 
modes that are not computationally or IO intensive, such as the 
incoherent station addition mode, are especially suited for paral- 
lel observing. This incoherent mode can therefore run commen- 
sally with many of the imaging projects (Rottgering et al. 2006; 
iFender eFa l. 2008) to continuously scan a large p art of the sky 
for in termittent pulsars such as PSR B 193 1+24 ( Kra mer et al. 
2006 ) and millisec ond radio transients (lLorimer et alj l200*7t 



Hessels et al. 2009). 



4.4. Survey strategies 

Three types of surveys were presented and their overall effi- 
ciency (Table [TJ is comparable. When deciding which LOFAR 
survey to implement, considerations relating to the different inte- 
gration times, beam sizes and operational requirements in these 
surveys may thus become more important: Core Coherent has 
short (order ~ 1 min) integrations, beneficial for finding acceler- 
ated systems and for system stability, and provides good local- 
isation. Superstation Coherent and Full Incoherent use of order 
~ 1 hr integrations, and are well suited for finding intermittent 
pulsars, but with only reasonable and poor initial localisation re- 
spectively. Core Coherent and Superstation Coherent produce 
much higher data rates than Full Incoherent. Tests on the differ- 
ent types of LOFAR data that are currently becoming available 
may potentially further differentiate between these strategies. 

5. Conclusions 

Because of its large area and field of view, LOFAR can reveal 
the local population of pulsars to a very deep luminosity limit. A 
25-day all-sky survey at 140 MHz would find 900 new pulsars, 
disclosing the local low-luminosity population and roughly dou- 
bling the number of pulsars known in the northern hemisphere. 
Millisecond pulsars in nearby globular clusters can be detected 
to lower flux limits than previously possible. Assuming the pul- 
sar population in other galaxies is similar to that in ours, we can 
detect periodicities or giant pulses from extragalactic pulsars up 
to several Mpcs away. 
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